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ABSTRACT 

The simultaneous presence of aerosol particles and condensible vapors In 

a saturated boundary layer may affect deposition rates to subcooled surfaces 

because of vapor-particle Interactions. Scavenging of condensible vapors by 

aerosol particles may lead to Increased particle size (and mass) and decreased 

vapor mass fraction, thus altering both vapor and particle deposition rates. 

Particles, If "sufficiently" concentrated, may also coagulate. Criteria are 

provided here to assess the significance of such phenomena when particles are 

already present In the mainstream and are not created Inside the boundary 

layer via homogeneous nucleatlon. It Is determined that there Is direct 

proportionality with (a) the mass concentration, (po>), of both condensible 

vapors and aerosol particles, and (b) the square of the boundary layer 

2 

thickness to particle diameter ratio, (6/d p ) . Inverse proportionality 

was found for mainstream to surface temperature difference, AT, If 

thermophoresis dominates particle transport. A numerical example Is worked 

out for a coal fired PFBC turbine application obtained from the literature. 

2 

It Is concluded that (S/d p ) Is the single most critical factor to 
consider In deciding when to neglect vapor-particle Interactions. 


*Case Western Reserve University, Department of Mechanical and Aerospace 
Engineering, Cleveland, Ohio 44106 and NASA Resident Research Associate. 
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NOMENCLATURE 


c p specific heat at constant pressure per unit mass of mixture 

D Brownian (Flck) diffusion coefficient (mass dlffuslvlty) 

Oa Oamkohler number, Eq. (12) 

dp particle mean diameter 

j" dlffuslonal mass flux 

k Boltzmann constant 

Le Lewis number, Pr/Sc 

Ma Mach number 

Nu Nusselt number for heat or mass transfer 

Pr Prandtl number 

p pressure 

R nose radius of turbine blade 

Re R Reynolds number based on blade nose radius 

r'“ local rate of vapor consumption, Eq. (1) 

Sc Schmidt number 

T absolute temperature 

t characteristic time 

Uoo turbine Inlet ("approach") velocity 

u,v fluid velocity In x and y directions parallel and normal to wall, 
respectively 

v^. thermophoretlc velocity, Eq. (3) 

x,y distance along and normal to the surface 

a thermal (heat) dlffuslvlty of carrier gas mixture, x/(pCp) 

^ thermal diffusion factor, Eq. (3) 

B binary coagulation rate constant, Eq. (14) 

6 boundary layer thickness 

X thermal conductivity of mixture 
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V 


viscosity of mixture 
v kinematic viscosity (momentum dlffuslvlty) of mixture, y/p 
p density 

w mass fraction In prevailing mixture 

Subscripts 
h heat transfer 

1 condensible vapor species 

m mass transfer 

max maximum 

p particle 

r reaction, Eq. (11) 

T thermophoresis 

w at the surface (wall) 

oo upstream Infinity 

HI scellaneous 
A difference operator 

INTRODUCTION 

The performance deterioration of many different types of high temperature 
systems Is closely related to the diffusion rates of reactants, reaction 
Intermediates, reaction products, and Impurities. Energy conversion device 
manufacturing as applied to aeronautics and/or power plants, and chemical 
synthesis reactors as applied to high technology materials processing (e.g., 
CVD techniques for microelectronic circuitry) are examples of such systems 
where deposition (via diffusion) Is suppressed or exploited to avoid 
efficiency degradation or to Improve product quality and life. 

Often both aerosol particles and condensible vapors experiencing steep 
temperature gradients are present In these systems. Here, we consider only 
sufficiently small particles that are Incapable of Inertlally Impacting on 
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deposition surfaces under the prevailing aerodynamic conditions. The recent 
studies of Ahluwalla and Im (1,2), Nadkarnl and Mahallngham (3), and Castillo 
and Rosner ( 4 , 5_) deal with the effects of having the simultaneous presence of 
vapors and particles and allowing for condensation In (super)saturated 
mainstreams and/or boundary layers. It has been shown that the mutual 
Interaction of vapors and particles has Interesting total (and Individual) 
deposition rate consequences and can significantly alter predicted behavior 
obtained from "pure" treatments alone. 

Small particle transport In temperature gradients Is dominated by 
thermophoresis (particle drift down a local temperature gradient) (6-9). 
Vapors, on the other hand, are transported primarily by Flck diffusion down 
local concentration gradients. Although under the prevailing thermodynamic 
conditions the supersaturation of the condensible vapors may not be sufficient 
to overcome the kinetic barriers for homogeneous nucleatlon, when both 
particles and vapors coexist In as system, particles can serve as possible 
nucleatlon sites for vapors. Vapors can be scavenged by condensation on the 
surface of particles. Moreover, If the particle concentrations and/or 
residence times Inside the boundary layer are large enough, the particles may 
coagulate (agglomerate). 

The objective of this study Is to provide a fast estimation method for a 
design engineer to assess the significance of (a) the Influence of the 
presence of condensible vapors on particle deposition rates, (b) the Influence 
of the presence of particles on vapor deposition rates, and (c) the Influence 
of coagulation on particle deposition rates, under conditions where the 
particles are already present In the mainstream and are not created Inside the 
boundary layer via homogeneous nucleatlon. An example case Is provided for a 
typical coal-fired pressurlzed-fluldlzed bed combustor (PFBC) power plant as 
obtained from Refs. 10 and 11. Note that the Intention Is to be able to 
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quickly decide when to safely Ignore the effects of above-mentioned 
phenomena. The reader Is referred to Refs. 1 to 5 If such effects can not be 
neglected. 

ANALYSIS 

The boundary layer conservation equations for condensible vapor species 
and small aerosol particles, even with such simplifying assumptions as 
Incompressible flow, constant (temperature Independent) thermodynamic and 
transport properties, very dilute species In the host gas (<V «p « 1). 
and monodlsperse aerosols, are coupled and have to be solved simultaneously. 
Thermophoresis dominates diffusion because even when the deposition surface Is 
about 2 percent cooler than the mainstream (T /T = 0.98 because of, 

W 00 

say, modest conduction and radiation losses with no active cooling), the 
enhancement of the deposition rate of, say, a 0.1 pm particle due to 
thermophoresis Is about an order of magnitude (even higher for larger 
particles and colder walls) (1). The boundary layer mass conservation 
equations for the condensible vapor and aerosol particles In two dimensions 
can then be expressed as: 
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where v^ Is the aerosol particle thermophoretlc velocity given by 


( 2 ) 
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(3) 


and 


1 II i 


denotes the mass of vapor locally scavenged by the aerosol particles 


per unit time and volume. If one approximates the shapes of the aerosol 
particles by spheres, neglects the relative velocity of condensible vapors 
with respect to aerosol particles, and takes the largest possible mass 
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transfer driving force, then (r"') ' can be expressed In terms of the 

diffusion flux of the condensible vapors onto the surfaces of aerosol 
particles, l.e., 



(4) 


where the mass transfer Nusselt number Is taken to be two, corresponding to 

mass transfer to/from spheres In a quiescent medium. Note that the analysis 

below will provide a conservative estimate of the significance of the 

above-mentioned effects because we will use the largest possible reaction 

(condensation) rate, (f 1 ") . In reality the actual reaction 

(condensation) rate, r'", will be smaller because (a) the local mass fractions 

of both condensible vapors, , and aerosol particles, o p , Inside the 

boundary layer will be smaller than their corresponding mainstream values 

«. and w , and (b) the driving force for condensation onto the 

1 ,e° P,« 

aerosol particle Is determined by the difference between the condensible vapor 
local mass fraction, <^ , In the boundary layer and the mass fraction locally 
prevailing on the aerosol particle surface at the gas/condensate Interface 
(assumed to be much smaller than here). 

The Effects of The Presence of Condensible Vapors on Particle Deposition Rates 
The thermophoretlc velocity, v.p of aerosol particles depends on the 
thermophoretlc dlffuslvlty, p D p , as shown In Eq. (3). However, because 
the thermophoretlc diffusion factor, <Xj p , Is Inversely proportional to the 
Brownian diffusion coefficient, D p , the thermophoretlc velocity Is actually 
Independent of D p (and, hence, d p ) (8). Therefore, It Is known a priori 
that vapors condensing on the aerosol particle surfaces will not affect their 
thermophoretlc velocity. However, the deposition rate of the aerosol 
particles will still be affected because the depositing particles will be 
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larger (heavier). A proper way to measure the significance of this effect Is 
to consider the maximum fractional change In the aerosol particle diameter due 
to vapor condensation. The rate of change of mass of an aerosol particle, 
Am/At, Is given by 


Am 

At = P p 



V 


Ad^ 
2 At 


(5) 


But Am/At Is also given by the mass deposition rate of the saturated vapors 
condensing on the aerosol particle surface, l.e.. 


Am 
At = 



Equating Eqs. (5) and (6), and solving for 


d p 

Ad 



gives 


( 6 ) 


Ad = 4 -4 — At 

d p 

P p P 


(7) 


In Eq. (7), At Is the residence time of the aerosol particles travelling with 
a thermophoretlc speed | v^. | Inside the thermal boundary layer of 
thickness, 6^, and can be estimated from 


It *h \ 

" V T ~ « T D * 

T.P P T w 4„ 


< 1.P U, P> 


w 

AT 


( 8 ) 


Combining Eqs. (7) and (8), the fractional change In the aerosol particle 
diameter Is given by 


Ad 

d 


Le 


1 


w *“1,- / w h 

• — 3 • — • I — — 


p ^T.p 16 ?^ r p \ p / 

The parameters In Eq. (9) can be readily estimated for the prevailing system 
operating conditions to assess the significance of vapor scavenging by 
particles (see the example below). 


AT 


(9) 
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The Effect of The Presence of Aerosol Particles on Vapor Deposition Rates 


A measure of the Importance of the presence of ash particles acting as 
sink for vapor molecules can be obtained by comparing the relative magnitude 
of the reaction ( r " 1 ) term In Eq. (1) to the diffusion term. An effective 
Damkohler number defined as 

_ _ characteristic vapor diffusion time across mass transfer boundary layer 

characteristic reaction (condensation) time 

Is such a quantity. The characteristic vapor diffusion time Is given by 



(10) 


and, the characteristic reaction (condensation) time Is given by 

i d p • 

r 12 D 1 • 

' 'max y ’ 

Therefore, 


(ID 


r P \ P , 


( 12 ) 


The quantities In Eq. (12) should be estimated for the prevailing system 
operating conditions. 

The Effect of Particle Coagulation on Aerosol Particle Deposition Rates 
Because we assume that aerosol particle transport Is dominated by 
thermophoresis, the residence time of the aerosol particle Inside the thermal 
boundary layer Is determined by Its thermophoretlc velocity. A comparison of 
the characteristic particle residence time Inside the thermal boundary layer 
to the characteristic coagulation time Is a measure of the significance of 
particle coagulation. If one assumes a nearly monodlsperse system with 
Brownian coagulation, then the ratio of two characteristic times Is given by, 
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2 


^thermophores Is 
^coagulation 




( 13 ) 


where 0 is the binary coagulation rate constant given by (12) . 


0 = 


8 kT 
3 v 


(14) 


Again, one can readily obtain the quantities required for Eq. (13) for the 
specific system in question (see also the example below). 


AN EXAMPLE CASE 


The example case in this study is taken from Refs. 10 and 11 for the 
Exxon PF8C pilot plant/cascade experiments. The plant operated at 10 atm with 
a combustion turbine inlet temperature of 871 °C (1144 K). The average 
particulate mass loading was 58 ppm and the particles had a density of 

3 

1.5 g/cm . The cumulative particle size distribution Is shown In Fig. 1 as 
duplicated from Ref. 10. The mean diameter of the particles can be read from 
Fig. 1 as about 1.2 pro, corresponding to the cumulative size distribution of 

3 

50 percent. Note that this mean diameter, based on particle mass (~d p ) 

rather than the more appropriate particle surface area (Sauter mean diameter), 

Is assumed to be sufficiently accurate for all practical purposes of this 

example. When Fig. 3 in Ref. 8 Is used as a guide to obtain the ratio of 

particle thermophoretlc dlffuslvlty to carrier gas thermal (heat) dlffuslvlty, 

a_ Le Is estimated to be about 0.2 for the prevailing system 
T,p P 

conditions. The maximum condensible vapor mass loading will be taken to be 
around 10 ppm and the condensible vapor will be reasonably assumed to have a 
Schmidt number of unity In the host combustion product carrier gas which Is 
basically air. The nose radius of the turbine blades used In the 
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plant was 0.635 cm. We will assume that the surface (wall) temperature of the 
blade is about 50 K (4 percent) smaller than the turbine Inlet temperature due 

to such losses as radiation, conduction, etc. If one Is Interested In the 

deposition rates at the stagnation point region of the blades, the heat and 
mass transfer Nusselt numbers can be estimated from Ref. 13: 

Nu h = J- = 0.81 (Re R ) 1/2 (Pr) 0,4 = 0.70 ReJ /2 (15) 

Nu m = = 0.81 (Re R ) 1/2 (Sc) 0 ' 4 = 0.81 ReJ /2 (16) 

m 

respectively, where the Prandtl number for air Is taken to be 0.7. Because It 
Is known that current turbines operate at subsonic velocities, we will assume 
that the turbine Inlet velocity Is around 250 m/sec resulting In a Mach number 
of about 0.4. We can now evaluate Eqs. (9), (12), and (13). The numerical 
values of each of the required quantities as well as their source references 
are given In Table I. The calculations show that 



Da = 5.8x10 

* thermophoresis _ ^ 8x1 0 _ ® 

^coagulation 

for the example case we have chosen. These values Indicate that all three 
effects considered here can be neglected for this example. 

The validity of the assumption that only convectl ve-dlf fuslonal mass 
transport Is Important and that the Inertial effects are absent can now be 
checked. The dimensionless Stokes number defined as Stk = (p d 2 U ) /( 1 8 yR) 

p p 00 

Is a measure of the effect of particle Inertia. The calculated Stokes number 
for this example Is 0.103, as compared with the critical Stokes number of 
0.125 for cylindrical target geometries above which aerosol particle 
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deposition (capture) starts to get dominated by Inertial Impaction. Note that 
although the example Stokes number Is less than the critical Stokes number, It 
Is sufficiently close to make It risky to completely neglect Inertial 
effects. However, Inclusion of such effects will result In Increased particle 
speeds, therefore reduced residence times, across the boundary layer thus 
minimizing the chance of Interaction of particles with condensible vapors or 
with each other. Hence, our conclusion that the simultaneous presence of 
aerosol particles and condensible vapors In the boundary layer does not affect 
the total (and/or Individual) deposition rate Is strengthened further. 
CONCLUSIONS 

When both aerosol particles and condensible vapors are simultaneously 
present In a saturated boundary layer, their Interaction may affect the 
Individual and/or total deposition rates to subcooled surfaces. Criteria are 
provided for neglecting the effect of aerosol particle scavenging of 
condensible vapors on (a) particle deposition rates, and (b) vapor deposition 
rates; and for neglecting particle coagulation. A numerical example Is solved 
for a typical coal fired PFBC turbine application obtained from the 
literature. It has been determined on the basis of numerical example that 
particle coagulation will hardly be a problem for turbine applications. 

Examination of Eqs. (9) and (12) reveals that the effect of particle 
scavenging of condensible vapors depends directly on the mass concentration of 
condensible vapors (pu, ) and aerosol particles (p« ), respectively. 

1,oo p,«o 

However, given the fact that the example was for p = 10 atm and that the mass 
fractions were high. It Is hard to conceive a real turbine application where 
the (pu) products will be an order of magnitude larger than those used In this 
example. The conservative estimation routine presented here, therefore, 
allows us to conclude that mass concentrations alone can not be responsible 
for vapor-particle Interactions for turbine applications. 
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Because most current (and future) turbines are actively cooled (larger 
AT's), the residence times of the aerosol particles Inside the boundary layer 
are. In fact, much smaller than the example case. This further reduces the 
vapor-particle Interaction effect on deposition rates (see Eq. (9)). 

A closer look at Eqs. (9), (12), and (13) shows that the ratio of the 

boundary layer thickness (note that the momentum, heat, and vapor species mass 

boundary layer thicknesses are essentially the same for our purposes) to the 

particle diameter Is the most Important factor In determining the significance 

of the scavenging and coagulation phenomena. There may certainly be some 

2 

cases In turbine applications where ( 6/d p ) can be significantly larger 
than the example case. This would result In nonnegllglble Ad p /d p 
and/or Da values. Lower mass flowrates (pU^) and trailing edge regions 
of larger size turbine blades may result In such Intolerably thick boundary 
layers. Slmllarlly, particle diameters may very well be smaller than, say, 

0.1 ym leading to large vapor condensation rates on existing aerosol particle 
surfaces. This observation Is especially Important If the boundary layer 
supersaturation Is so large that homogeneous nucleatlon (not treated here), 
may Initiate, because the resulting particle diameters from homogeneous 
nucleatlon are usually between 0.01 to 0.1 ym, leading to significant 
vapor-particle couplings. The reader Is referred to Refs. 1 to 5 for the 
treatment of such coupled transport phenomena. 

Finally, one should note that the analysis presented here applies also to 
turbulent boundary layers as long as thermophoresis Is the dominant particle 
transport mechanism or the analysis can easily be extended to turbulent 
boundary layers using the same concepts but substituting turbulent momentum, 
heat, and mass dlffuslvltles. 
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PARTICLE DIAMETER, |im 


TABLE I. 


Symbol 

Value 

Source 

c p 

0.2787 cal/gr-K 

Ref. 14 

dp 

1.2 vm 

Ref. 10 

k 

1 .38x10-23 J/K 

Ref. 14 

Le* = Pr/SCA 

0.7 

Calculated 

Ma 

0.4 

Assumed 

Nu h 

229 

Calculated 

Nu m 

265 

Calculated 

P 

10 atm 

Ref. 10 

Pr 

0.7 

Ref. 14 

R 

0.635 cm 

Ref. 10 

Rcr = (pUcoR)/y 

1 .07x10® 

Calculated 

Sc* 

1 

Assumed 

To, 

1144 K 

Ref. 10 

AT 

50 K 

Assumed 

T W 

1094 K 

Calculated 

Uca 

250 m/sec 

Calculated 

a = \/(pC p ) 

0.21 cm 2 /sec 

Calculated 

a T.pL e p 

0.2 

Ref. 8 

G = (8^/(3 y) 

9.2x10-11 cm 3 /sec 

Calculated 

5 h 

27.7 ym 

Calculated 

5 m 

24.0 ym 

Calculated 

\ 

179xl0" 6 cal/cm-sec-K 

Ref. 14 

y 

460xl0 -6 gr/cm-sec 

Ref. 14 

p 

3.1x1 0 -3 gr/cm 3 

Ideal gas law 

p P 

1 .5 gr/cm 3 

Ref. 10 


10 ppm 

Assumed 

Wp,oo 

58 ppm 

Ref. 10 



CUMULATIVE PERCENT GREATER THAN 


Fig. 1. - Size distribution for 10:1 pressure ratio turbine. 
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of the Automotive Stirling Engine 
Joseph R. Stephens 
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Abstract 

The automotive Stirling engine now under development by DOE/NASA as an 
alternative to the Internal combustion engine, Imposes severe materials 
requirements for the hot portion of the engine. Materials selected must be 
low cost and contain a minimum of strategic elements so that availability Is 
not a problem. Heater head tubes contain high pressure hydrogen on the Inside 
and are exposed to hot combustion gases on the outside surface. The cylinders 
and regenerator housings must be readily castable Into complex shapes having 
varying wall thicknesses and be amenable to brazing and welding operations. 
Also, high strength, oxidation resistance, resistance to hydrogen permeation, 
cyclic operation, and long-life are required. A research program conducted by 
NASA Lewis focused on alloy chemistry and mlcrostructural control to achieve 
the desired properties over the life of the engine. Results of alloy 
selection, characterization, evaluation, and actual engine testing of selected 
materials will be presented. 

KEY WORDS: Automotive Stirling engine; Alloy chemistry; Microstructure; Creep 

resistance; Oxidation; Hydrogen permeation; Iron-base alloys 



The Department of Energy (DOE) and NASA [1] are currently Investigating 
the Stirling engine as an alternative power source to the Internal combustion 
engine for automotive applications. The Stirling engine (Fig. 1) Is an 
externally heated engine that offers the advantages of high efficiency, low 
pollution, low noise, and multifuel capability. Hydrogen Is used as the 
working fluid for automotive applications In order to obtain maximum 
efficiency from the engine. Heat Is Input to the working fluid by a 
combustion flame Impinging on the outside walls of the tubes containing the 
hydrogen working fluid. Current engines designed for automotive applications 
use tubes with a wall thickness of 0.75 mm (3.0 mm l.d. by 4.5 mm o.d.). The 
thin wall Is required to achieve maximum heat Input to the hydrogen. However, 
the tube material must be of sufficient strength to contain the hydrogen, 
which operates at high pressures and temperatures. Automotive design 
requirements call for a normal maximum pressure of 15 MPa, with peak pressures 
reaching 21 MPa during transient conditions. Tube walls are designed for an 
average temperature of 820 °C, with 870 °C common on the flame side. A 
further requirement Is that the tubes operate for a period of 3500 h, which Is 
approximately equivalent to the 160 000-km (100 000-mlle) driving life of an 
engine. 

The tubes currently are heated by direct exposure to a combustion flame 
fired by gasoline or diesel fuel. Alternative fuels such as alcohol and shale 
oil derivatives may be considered for future applications. Oxidation and 
corrosion resistance Is one of several criteria for selecting heater-head-tube 
alloys to withstand the extreme operating requirements of the Stirling 
engine. A second criterion Is the resistance of the tube material to hydrogen 
permeation so that hydrogen recharging of the engine Is required In excess of 
6 month Intervals. Finally, the creep-rupture strength for 3500 h life Is a 
determining factor In selecting a heater tube alloy. These criteria are 
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Interrelated. For example, excessive scale spalling due to 
oxidation/corrosion will reduce the tube wall thickness and lead to premature 
failure caused by the high-pressure hydrogen rupturing the tubes or permeating 
rapidly through the tube walls. Prototype engines and Initial models of 
experimental automotive engines use N-155 (Multlmet) for heater head tubes. 
However, this Iron-base alloy contains 20 percent cobalt, a costly and 
strategic metal. Efforts are currently under way by NASA [2-7] and the 
contractors on this program, Mechanical Technology Inc. (MTI) and United 
Stirling Sweden (USAB), to Identify substitutes for N-155. The primary 
emphasis Is on high-strength austenitic Iron-base alloys although a limited 
number of nickel-base superalloys are under consideration. As part of the 
evaluation of candidate substitute alloys the oxidation and corrosion 
resistance has been determined for 16 alloys [8]. Hydrogen permeation 
coefficients have been determined for most of those alloys [9], and strength 
properties have been measured by creep testing In air [2] and hydrogen [10]. 
Further, endurance testing In a materials simulator rig and engine testing 
have been used to characterize selected candidate heater tube alloys [3,6,9]. 

The purpose of this paper Is to describe the effects of alloy chemistry 
and microstructure on the oxidation/corrosion behavior and resistance to 
hydrogen permeability of candidate alloys. 

Procedure 

Materials 

Twenty alloys (14 Iron base, five nickel base, and one cobalt base) were 
evaluated for their oxidation and corrosion resistance, hydrogen permeability, 
and endurance under conditions of temperature and environment anticipated for 
automotive Stirling engines. The chemical compositions of the alloys are 
listed In Table 1. All of the alloys are commercially available except 
AL-EX-20, an experimental alloy supplied by Allegheny Ludlum Steel Company. 
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Iron-base alloys are the primary candidates for automotive Stirling engine 
applications because their raw material costs are lower than those for 
nickel-base and cobalt-base alloys. Nickel-base alloys are considered as 
possible candidates In the event a less expensive Iron-base alloy cannot be 
Identified that will meet all of the requirements of the Stirling engine. 
Cobalt-base alloys are not under consideration for automotive Stirling engine 
applications because of the high cost of cobalt. Its limited availability, and 
Its status as a strategic material (more than 95 percent of the cobalt used In 
the United States Is Imported). However, one cobalt-base alloy (HS-188) was 
Included In this Investigation as a comparison with the Iron- and nickel-base 
alloys . 

For oxidation testing, test coupons 1.27 by 2.54 cm, and generally 0.8 to 
1.6 mm thick, with a 0.3-cm-dlameter hanger hole, were cut from the alloy 
sheet material. For hydrogen permeability and endurance testing the alloys 
used In this study were obtained commercially In the form of tubing with an 
outside diameter of 4.5 mm and an Inside diameter of 3.0 mm (the size used In 
current prototype engines). Four of the tubing alloys were weld-drawn and the 
others were seamless tubing. The weld-drawn alloys were N-155, 19-9DL, 

Inconel 718(wd) and HS-188. 

Stirling Engine Simulator Materials Test Rig 

The Stirling engine simulator rig used In this program was designed and 
fabricated at the NASA Lewis Research Center; It consists primarily of a 
combustion gas heating chamber (Fig. 2) with auxiliary heating, control, and 
gas management systems. The rig and Its operation have been described In 
detail [3]. The rig was fired with natural gas for Initial Ignition and then 
with diesel fuel throughout the test run. Duplicate oxidation-corrosion test 
specimens were suspended by platinum wires from the hairpin turns of the tubes 
shown In Fig. 2. The specimens were heated to the desired 820 °C test 
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temperature (automotive Stirling engine temperature) by the combustion flame 
flowing over them. The rig was used primarily as a means of ranking candidate 
heater-head- tube alloys [3,6,9]. The oxidation-corrosion testing was 
conducted simultaneously with a 3500-h endurance and permeability run of 
tubing alloys. An endurance run consisted of a series of 5-h cycles to obtain 
the required 3500 h of test time to simulate the minimum life of an automotive 
engine. A typical heating cycle consisted of a 6- to 10-mln preheat to get to 
the operating temperature, a 5-h hold at temperature, and a 1-h cooldown to 
near room temperature. The change In specific weight AW/A with time t for 
the oxidation coupons was fitted by least squares [11] to equation (1). 


f . k ; /2 t i/2 . * 2 t 1SEE 


Hi 


where k^ Is an oxide growth constant comparable to a parabolic oxidation 
scaling constant, k,, Is an oxide spalling constant, and SEE Is the 
standard error of estimate. An oxidation attack parameter K [11] was 
derived from equation (1) and Is defined as 


1/2 

K a = k.j * 10k 2 (2) 

The original hydrogen pressure In the tubes and the pressure at the end of a 
5-h cycle were used to calculate the hydrogen permeability coefficient, <t>, 
using the equation: 


where 



<|>AP s Tt 

2LVT 

s 


(3) 


P pressure In closed system, MPa 

P original pressure, MPa 

o 

A permeated area, cnf* 

P $ standard pressure, MPa 

T temperature of system, K 
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t time, sec 

L membrane thickness, cm 

3 

V volume of system, cm 

standard temperature, K 

Results 

Stirling Engine Simulator Oxidation and Corrosion Tests 

Weight Change. - The alloys could be categorized Into four groups on the 

basis of specific weight change data. Alloys In group I had excellent 

oxidation and corrosion resistance, and their specific weight change was 

positive (weight gain) throughout the 3500-h test. Group I comprised the five 

alloys CG-27, Incoloy 800, HS-188, N-155, and Inconel 718. Three of the 

alloys, CG-27, Incoloy 800, and HS-188, did not exhibit any significant 

spalling during testing In the Stirling engine simulator rig. Final average 

2 

specific weight change (Fig. 3) ranged from 1.4 to 0.17 mg/cm for CG-27 to 
Inconel 718, respectively. 

A second group of alloys Is considered to have relatively good oxidation 
and corrosion resistance under simulated Stirling engine operating 
conditions. The five alloys In group II In order of Increasing final weight 
loss are RA-330, Sanlcro 31H, Sanlcro 32, 12RN72, and 253 MA. Final average 

2 

specific weight change (Fig. 4) ranged from -0.54 to -6.0 mg/cm for RA-330 
to 253 MA, respectively. 

The third group of alloys, 316 stainless steel, W-545, A-286, and 

Nltronlc 40 Is considered to have fair to poor oxidation and corrosion 

resistance on the basis of specif Ic-welght-change-versus-tlme data. Final 

2 

average specific weight change (Fig. 5) ranged from -22.8 to -41.8 mg/cm 
for 316 stainless steel to Nltronlc 40, respectively. 

A final group of alloys, consisting of 19-9DL and AL-EX-20, exhibited 
catastrophic oxidation and corrosion behavior upon rig testing. Spalling 
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began after only a very short time Into the test and excessive weight loss 
occurred after 500 h of cyclic testing. The tests were terminated before the 
3500-h goal because of the unacceptable oxidation and corrosion resistance of 
these alloys. 

Attack parameter K fl . - Oxidation and corrosion 
welght-change-versus-tlme data were fitted to the parallnear equation 
(eq. (1)). For three of the alloys, HS-188, Incoloy 800, and CG-27, where no 
spalling was observed, the attack parameter was calculated from equation (2) 
with k 2 = 0. The remaining 13 alloys also derived their attack parameters 
from equation (2) but with values of k^ Included. For alloys 19-9DL, 
Nltronlc 40, and AL-EX-20 the k^ 2 term dropped out so that 
K = 20k o t was used. Figure 6 Is a summary plot of calculated K 
values for the alloys. Values of K are plotted on a log scale In 

3 

decreasing order of oxidation and corrosion resistance based on the highest 
values of the duplicate specimens. The attack parameter ranged from 0.011 for 
HS-188 to 12.49 for 19-9DL, or over three orders of magnitude. Examination of 
the K values In Fig. 6 suggested that the alloys could again be 

3 

categorized Into four groups. Group I comprised the alloys HS-188, Sanlcro 
31 H, Inconel 718, Incoloy 800, N-155, CG-27, and RA-330, with K ranging 

3 

from 0.011 for HS-188 to 0.031 for RA-330. Group II comprised the three 
alloys Sanlcro 32, 12RN72, and 253 MA, with K values of 0.098, 0.110, and 

3 

0.159, respectively. Group III comprised the alloys A-286, W-545, 316 
stainless steel, and Nltronlc 40, with K ranging from 0.370 to 0.659. 

3 

Group IV comprised AL-EX-20,, with a K of 4.68, and 19-9DL, with a K 

3 3 

of 12.49. This grouping was somewhat arbitrary, but corresponded to the 
previous grouping based on final weight change. The only exception was for 
RA-330, which experienced a net weight loss and was placed In group II 
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(Fig. 4) but has a K value similar to those for the other five alloys In 

a 

group I. 

X-ray diffraction. - The x-ray diffraction data from the specimen surface 
after testing are summarized In Table 2 by rank of the alloy In decreasing 
X-ray Intensity. In general, the most oxidation and corrosion resistant 
alloys formed nearly pure Cr^O^ as either the most or secondmost Intense 
oxide. The Cr 2 0 3 oxide Is listed as a sesquloxlde, d (024) = 0.1818 nm In 
Table 2. Because Fe went Into solution as Fe 2 0 3 , the value of d 
Increased. As the d value approached 0.1839 nm, the oxide became 
essentially pure Fe 0 . Chromite spinel with a ranging from 0.835 to 

L J 0 

0.840 nm also formed on most of the alloys during oxidation and corrosion 
testing. This chromite spinel appeared to be Innocuous. Alloy A-286 showed 
another spinel with a Q = 0.830 nm along with a lower sesquloxlde with d = 
0.1813 nm, and also N10 was formed. Alloy 19-9DL, which oxidized 
catastrophically, formed as Its two strongest oxides a spinel with a Q = 

0.840 nm and the desirable oxide Cr 2 0 3 with d = 0.1821 nm. 

Metallography. - Microstructures of selected alloy specimens after the 
3500-h test are shown In Fig. 7. Alloy C6-27 (Fig. 7(a)) had an adherent 
oxide scale on the specimen surface. Alloy Inconel 718 (Fig. 7(b)) had a 
thin, adherent oxide scale with a small depletion zone present In this alloy. 
Sanlcro 32 (Fig. 7(c)) formed large voids In the depletion zone. Severe grain 
boundary attack occurred In W-545 (Fig. 7(d)) and resulted In almost total 
consumption of the alloy matrix beneath the surface oxide. 

Electron microprobe analysis. - Alloy CG-27, a group I alloy, was 
examined because It had excellent oxidation resistance and contained the least 
chromium of the 16 alloys tested In the rig. The distribution of the major 
alloying elements and oxygen In the oxide scale and the adjacent metal 
substrate of CG-27 (Fig. 8) showed the surface oxide to be rich In Iron, 
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nickel, and to a lesser extent aluminum. Beneath this oxide layer was a thin, 
continuous oxide layer rich In chromium and titanium. Adjacent to the 
chromium-rich Intermediate scale was an aluminum-rich oxide next to the metal 
substrate, which exhibited a wide zone of Internally oxidized aluminum-rich, 
and somewhat fewer titanium-rich, protrusions. 

Metal lographlc results for the group III alloy W-545 (Fig. 7(d)) 

Indicated a large surface scale consisting of a porous oxide. Results from 
electron microprobe analyses (Fig. 9) showed an outer oxide rich In Iron, 
nickel, and chromium. Oxygen had penetrated uniformly Into this alloy to a 
substantial depth. The oxide formed appeared to be rich In nickel with some 
chromium. Iron was depleted from this area, but discrete areas of molybdenum 
were observed near the surface. Aluminum and titanium did not appear to enter 
Into the oxidation process of this alloy. The composition of W-545 Is similar 
to that of CG-27 except for the large difference In aluminum content. This 
suggests that aluminum makes a major contribution to the oxidation resistance 
of CG-27 and that Its low concentration In W-545 may explain the poor 
oxidation resistance of that alloy. 

Stirling Engine Simulator Hydrogen Permeation and Endurance Tests 

Hydrogen permeability. - Previous measurements of hydrogen permeation 
during cyclic testing In the simulator rig [12] Indicated that hydrogen loss 
was rapid through the thin wall tubes of alloy N-155. Hydrogen pressure decay 
curves are shown In Fig. 10(a) for selected 5 h cycles In an 100 h test. It 
should be noted that pressure loss Is rapid and Is most severe for the first 
5 h cycle. The slight Improvement for the remaining cycles Is attributed to 
the formation of an oxide scale on the exterior surface of the tubes. It was 
shown previously [13] that a thin (0.02 mm) carbide could form on the tubes 
Internal surface In an engine presumably due to oil leaking Into the heater 
head. This carbide layer resulted In reducing the permeability rate of 
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hydrogen from the engine. A program was undertaken at NASA Lewis [12] to 
Investigate various dopants such as CO, CO,,. CH^, etc. to form carbide or 
oxide surface layers on the tubes Internal surfaces to act as a barrier to 
hydrogen permeability. The effects of 1 vol % CO^ In hydrogen are shown 
In Fig. 10(b) for alloy N-155. The results suggest that the oxide that formed 
on the Internal surfaces was effective in reducing hydrogen loss from the 
tubes. Hydrogen permeability coefficients were calculated from the pressure 
decay curves using equation (3). Permeability coefficients as a function of 
rig exposure time are plotted In Figs. 11 to 13 for Iron, nickel, and 
cobalt-base alloys, respectively. For Iron-base alloys, hydrogen permeability 
coefficients vary over two orders of magnitude while for the five nickel-base 
alloys a variation of about a factor of two Is noted. A comparison of the 
hydrogen permeability coefficients after 250 h testing (since some alloys did 
not survive much longer) Is shown In Fig. 14. In general the nickel-base 
alloys exhibited lower permeability coefficients than the Iron-base alloys. 

An Iron-base alloy, CG-27 was an exception to this general observation. 

Endurance testing. - Endurance testing was also conducted on the tube 
alloys by filling with 15 MPa helium. Since helium did not permeate the tube 
walls, a constant pressure (stress) was maintained. This test Is much more 
severe than actual engine operation where the mean pressure Is near 4 MPa, but 
does provide an effective method for ranking the strengths of the alloys under 
simulated engine conditions. A comparison of rupture lives Is shown In 
Fig. 15. Only CG-27 and Pyromet 901 survived the 3500 h test. It should be 
noted for alloy Inconel 718 that seamless tubing (Inconel 718(a)) had a much 
longer life than the weld-drawn tubing (Inconel 718(wd)). Cold working of 
12RN72 (cw) reduced the rupture life from the standard heat treated condition. 

Metallography. - Microstructures of selected alloys are shown In 
Fig. 16. Alloy CG-27 Is shown before and after exposure. The oxide film on 
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both Internal and exterior surfaces Is evident for the hydrogen and 1 percent 
C0 2 containing tube (Fig. 16(a)). Alloy Inconel 625 (Fig. 16(b)) exhibited 
external scale cracking after cyclic exposure. 

Discussion of Results 

Some Insight Into the excellent oxidation and corrosion behavior of CG-27 
can be gained from the electron microprobe results. The chromium distribution 
shows that CG-27 had a chromium-rich oxide layer ranging In thickness from 
about 3 to 4 pm. The poorer oxidation behavior of group II and III alloys was 
characterized by chromium-rich oxide layers about 8 pm thick (group II alloys) 
and from 20 to 35 pm thick (group III alloys) [8]. All 16 alloys studied 
herein formed Cr 0 chromite spinels as their major oxides. The thin, 

C «J 

adherent oxide formed by the group I alloys CG-27 and Incoloy 800 (not shown) 
Is believed to be a major contributor to the oxidation resistance of these 
alloys . 

Further comparison of the alloys characterized by electron microprobe 
techniques Indicated that Iron and nickel were present as a thin outer oxide 
In group I alloys. In contrast. In group II alloys and especially In 
group III alloys. Iron and nickel were observed as thick, discrete areas of 
oxide that were In some cases ready to spall from the specimen. Evidently In 
group I alloys the Iron and nickel oxides formed early In the 3500-h test, and 
then oxide formation of the more reactive alloying elements took over the 
oxidizing process. The thicker oxides and higher spalling rates In particular 
for group III alloys (specific weight losses of 23 to 42 mg/cm as compared 
with little or no weight loss for group I alloys) suggest that Iron and nickel 
oxides formed throughout the 3500-h test. 

Particular Interest was focussed on the apparent roles of the more 
reactive oxide-forming alloying elements aluminum and titanium contained In 
some of the alloys. In CG-27 Internal oxidation of aluminum to a depth of 
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approximately 10 pm below the oxide-metal Interface was evident. In addition, 
aluminum was present In the surface oxide, and a thin oxide layer rich In 
aluminum existed adjacent to the metal substrate. A similar, though not as 
evident, distribution of aluminum occurred In Incoloy 800. Incoloy 800 
contains 0.38 wt % aluminum as compared with 1.5 wt % In CG-27. A thin, 
discontinuous aluminum-rich oxide layer near the oxide-metal Interface and 
grain boundary oxides were previously noted for Sanlcro 32. Titanium, at a 
concentration of 2.5 wt % In CG-27, appeared to be associated with aluminum 
In the oxide scales and as Internally oxidized. In contrast, for the 
group III alloy W-545, which contains slightly more titanium (2.85 wt %) but 
much less aluminum (0.2 wt %) , titanium does not preferentially oxidize. 

The role of reactive elements such a yttrium, scandium, cerium, 

zirconium, and hafnium In Improving oxide scale adherence particularly for 

Fe-Cr-AL, Nl-Cr-Al, and Co-Cr-Al alloys has been Investigated In much detail 

[14-19]. The standard free energies of oxide formation of these reactive 

elements are more negative than those for Iron, nickel, cobalt, chromium, and 

aluminum. Small additions (<1 at %) of these elements have been shown to 

dramatically Improve oxide scale adherence and to reduce spalling In cyclic 

oxidation. The exact mechanism by which such an Improvement occurs Is not 

agreed upon by Investigators In this field. For example, yttrium and scandium 

have been reported to Increase oxide growth rates by forming rapid diffusion 

paths along yttrla and scandla stringers [14]. Increased grain boundary 

18 

transport of oxygen has been demonstrated by 0 tracer studies on Nl-Cr-Al 
alloys with zirconium additions [15]. The resulting oxide adherence has been 
attributed [14,15,18] to (1) mechanical oxide pegging; (2) stress-strain 
accommodations at the oxide-matrix Interface; (3) Increased chemical bonding 
across the Interface; and (4) a vacancy sink mechanism that precludes void 
formation at the oxide-metal Interface, which would promote spalling. 
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It Is postulated that for the alloys studied herein, aluminum plays a 
role similar to that of the reactive elements discussed previously. Aluminum 
has a more negative standard free energy of oxide formation than do the major 
alloying constituents of the candidate heater-head-tube alloys such as Iron, 
nickel, and chromium. The Internal aluminum oxidation suggests that 
controlling the Inward oxygen diffusion plays a major role In the oxidation of 
CG-27 and probably Incoloy 800 [8]. Although the remaining group I alloys 
were not studied by electron microprobe techniques, the reactive elements 
niobium In N-155, niobium and aluminum In Inconel 718, and lanthanum In HS-188 
probably play a role similar to that of aluminum In the two alloys that were 
examined. The results further suggest that although 0.4 wt % aluminum Is 
adequate to be effective In promoting good scale adherence (Incoloy 800, 
Sanlcro 32, Sanlcro 31H, and Inconel 718), a concentration of 0.2 wt % 
aluminum such as exists for W-545, Fig. 9 and A-286 Is Insufficient to Improve 
the oxidation behavior of these alloys. It Is not clear from this study which 
of the previously suggested models for oxide scale adherence Is to be 
preferred . 

The factors effecting oxidation/corrosion behavior of the alloys Is also 
Important In minimizing hydrogen permeation through the tube walls. Formation 
of an adherent oxide film on the Internal surface of the heater head tubes Is 
required to reduce the hydrogen loss by permeation. Alloys with low 
oxidation/corrosion attack parameters, K also had low hydrogen 
permeability coefflcents, 4 > . Results presented herein have clearly shown 
the Importance of alloy composition, specifically the presence of the more 
reactive elements such as A1 , T1, and La as a key to controlling 
oxidation/corrosion behavior of the alloys. This, In turn. Is reflected In 
the resulting microstructure where thin oxide layers containing the reactive 
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elements are desirable for long-life, low-hydrogen permeability heater head 
tube alloys. 

Application of Results 

Based on the results of these tests alloys CG-27, 12RN72, Sanlcro 32, and 
Inconel 625 were selected for engine testing. Engine testing under 
accelerated conditions exceeded 2000 h or an equivalent 3300 h of 
urban-highway driving. Only alloy CG-27 did not fall by creep rupture In the 
engine test [20]. Alloy CG-27 has been selected as the heater head tube 
material for the automotive Stirling engine because of Its good creep-rupture 
(endurance) strength, resistance to oxidation and corrosion, and resistance to 
hydrogen permeation as a result of forming a thin, tenacious oxide using C0 2 
doped hydrogen. 

Concluding Remarks 

This study has shown that chemistry and microstructure of heater head 
tube alloys which develops during operation of the engine and In simulated 
tests play an Important role In alloy selection for the automotive Stirling 
engine. Specifically the following features were shown to be Important: 

1. The presence of a sufficient amount (about 0.4 wt %) of the reactive 
metal aluminum In Iron-base alloys leads to a microstructure characterized by 
a thin, adherrent oxide scale rich In aluminum plus Internal oxidation of 
aluminum which Impedes Inward diffusion of oxygen. 

2. The thin, adherrent oxide scales on both the combustion gas side and 
the CO^ doped hydrogen side of the tubes act as barriers to hydrogen 
permeation. 

3. The thin, adherrent oxides characterized by no or very little spalling 
and metal loss lead to utilization of the full cross-sectional area of the 
tube alloy to provide the required strength for the 3500 hr life of an engine. 
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TABLE 1 — CHEMICAL COMPOSITION OF COMMERCIAL ALLOYS 



^Evaluated in hydrogen permeability and endurance testing only. 
^Evaluated in oxidation/corrosion testing only. 





















TABLE 2 — X-RAY DIFFRACTION DATA OF SPECIMEN SURFACE AFTER 3500-h TEST 
IN ORDER OF DECREASING INTENSITY 
[Values of d and a Q are in nanometers.] 
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Figure 1. - Schematic representation of automotive Stirling engine, 



Figure 2. - Schematic representation of Stirling engine materials simulator test rig. 
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Figure 3. - Specific weight change data for 
group I alloys tested at 820 °C for 
3500 hr. 
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Figure 4. - Specific weight change data for 
group II alloys tested at 820 °C for 3500 hr. 
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Figure 5. - Specific weight change data for 
group III alloys tested at 820 °C for 
3500 hr. 
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Figure 7. - Photomicrographs of selected 3500 hr oxidation and corrosion test specimens. 
Longitudinal view showing surface scales and depletion zones; magnification, 250. 
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Figure 8. - Concluded. 
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Figure 9. - Concluded. 
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Figure 10. - Pressure decay curves for 
N-155 at 820 °C. 
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Figure 11. - Hydrogen permeability coefficient <t> versus rig exposure time for iron 
base alloys tested at 820 °C with H 2 + 1% CO 2 at 15 MPa. 
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Figure 12. - Hydrogen permeability coefficient 
<t> versus rig exposure time for nickel base 
alloys tested at 820 °C with H 2 + 1% C0 2 at 
15 MPa. 
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Figure 13. - Hydrogen 
permeability coefficient 
0 versus rig exposure 
time for cobalt-base 
alloy H S - 188 tested at 
820 °C with H 2 + 1% 
C0 2 at 15 MPa. 
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Figure 14. - Apparent hydrogen permeability coefficient at 250 hr for tubing material endurance tested at 820 °C and 15 MPa. 
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Figure 15. - Ranking of hairpin tubes according to rupture lives when pressurized with 15 MPa helium at 820 °C. 






Unexposed 


Unexposed 


H 2 + 1% C0 2 exposed 


H 2 - 1% C0 2 exposed 


(b) Inconel 625. 


Figure 16. - Microstructures of hairpin tubes before and after rig exposure, 820 C, 3500 h. 
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